New Insights into the Inhibition of Hesperetin on Polyphenol Oxidase: Inhibitory Kinetics, Binding Characteristics, Conformational Change and Computational Simulation

The inhibitory activity of hesperetin on polyphenol oxidase (PPO) and their interaction characteristics were investigated using multiple spectroscopic methods and computational simulation. Hesperetin, a mixed inhibitor, reversibly inhibited PPO activity, and its half-maximum inhibitory concentration (IC50) values on monophenolase and diphenolase were 80.8 ± 1.4 μM and 776.0 ± 15.5 μM, respectively. Multivariate curve resolution–alternate least squares (MCR–ALS) analysis suggested PPO interacted with hesperetin and formed PPO–hesperetin complex. Hesperetin statically quenched PPO’s endogenous fluorescence, and hydrophobic interactions mainly drove their binding. Hesperetin affected the polarity of the microenvironment around the Trp residues in PPO, but had no effect on that around Tyr residues. Circular dichroism (CD) results showed that hesperetin increased α-helix content and decreased β-fold and random coil contents, thus tightening PPO’s structure. Molecular docking showed that hesperetin entered the hydrophobic cavity of PPO, bound near the dinuclear copper active center, interacted with Val283, Phe264, His85, Asn260, Val248, and His263 via hydrophobic interactions, formed hydrogen bonds with Met280, His89, and His259 residues and also interacted with Phe292, His61, Phe90, Glu256, His244, Asn260, Phe264, and Gly281 via van der Waals forces. The molecular dynamics simulation results also demonstrated that the addition of hesperetin reduced the stability and hydrophobicity of PPO and increased PPO’s structural denseness. Thus, the inhibition of hesperetin on PPO may be because hesperetin bound near the active center of PPO, interacted with the surrounding residues, occupied the binding site for substrate, and induced the changes in PPO’s secondary structure, thus inhibiting the catalytic activity of PPO. This study may provide novel views for the inhibition of hesperetin on PPO and theoretical guidance for developing flavonoids as new and efficient PPO inhibitors.


Introduction
Polyphenol oxidase (PPO), also called tyrosinase, can cause the browning of vegetables, fruits and shrimps in the food industry [1]. As a rate-limiting enzyme related to enzymatic browning reaction, it catalyzes monophenols to diphenols and further oxidizes diphenols to o-quinone, finally transforming to melanin [2]. The conversion of phenols to melanin not only decreases the nutrition of foods and produces some harmful substances, but also causes food to show an unacceptable color and lose commodity value [3]. Thus, it is worthwhile to inhibit the activity of PPO to avoid food waste.
The traditional ways to reduce the activity of PPO in food can be divided into physical and chemical methods. Physical methods usually include heating treatment, high-pressure treatment and UV treatment, and chemical methods usually involve adding acidifiers (citric acid, ascorbic acid, etc.), or sulfur dioxide and its sulfites [4]. However, when the above methods inactivate PPO, they may also destroy the nutrients in food. For instance, PPO (from mushroom, EC 1.14.18.1, 128 kDa) was obtained from Worthington Biochemical Co. (Lakewood, CO, USA). Hesperetin (purity ≥ 98%), L-tyrosine, and L-dopa were obtained from Aladdin Biochemical Technology Co. (Shanghai, China). Kojic acid (KA, purity > 99%) was acquired from Oddfoni Biological Technology Co. (Nanjing, China). Other chemicals used were of analytical grade.

PPO Activity Assay
To determine the monophenolase inhibitory activity of hesperetin, PPO (0.2 µM) was mixed with various concentrations of hesperetin (0, 25,50,100,150,200, 250 and 500 µM), and the mixtures were incubated at 30 • C. After 1 h, L-tyrosine (0.5 mM, substrate) was added to the mixtures, and the detection began. The absorbance values were recorded at 475 nm for 30 min (121 times per 15 s) using a microplate reader (Varioskan LUX, Thermo Scientific, Waltham, MA, USA) to determine the oxidation of L-tyrosine [18].
Based on a previous method [19], the diphenolase inhibitory activity of hesperetin was determined spectrophotometrically by measuring the rate of L-dopa converted to dopaquinone. A series of assay solutions containing 20 µL PPO (20.0 µM) and various amounts of hesperetin (20.0 mM) were incubated at 30 • C. After 1 h, 200 µL L-dopa (5.0 mM, substrate) was added to start the action. The absorbance values were continuously monitored at 475 nm for 200 s (41 times per 5 s) by a UV-vis spectrophotometer (UV-2450, Shimadzu, Tokyo, Japan). Kojic acid was used as a positive control [20]. The inhibitory capacity of samples with different concentrations was expressed through relative PPO activity. The PPO activity without an inhibitor was defined as 100%, and the relative activity of PPO suppressed by different samples; relative activity (%) = K/K 0 × 100%, where K and K 0 represent the slopes of the absorbance, changed over time with and without an inhibitor, respectively [21].

Kinetic Analysis of Inhibitory Type
The inhibitory type of hesperetin against PPO could be evaluated using Lineweaver-Burk plots, and the mixed inhibition was described through Equation (1) [22]: where v and V max represent enzyme reaction velocity and its maximum value, respectively.
[I], [S] and K m are the concentrations of hesperetin and L-dopa, and the Michaelis-Menten constant, respectively. α denotes the apparent coefficient and K i represents dissociation constant of hesperetin binding with free PPO, which can be calculated using the secondary Equations (2) and (3) [23].

Multivariate Curve Resolution-Alternate Least Squares (MCR-ALS) Algorithm
MCR-ALS is acknowledged as a robust chemometric methodology to analyze massive and sophisticated data [24], and was used to obtain the concentration curves and spectral information about hesperetin, PPO and the PPO-hesperetin complex. Fluorescence titration experiment A: The concentration of PPO was 4.0 µM and the same volume of hesperetin was added into PPO solution 25 times (0→20.0 µM). Fluorescence titration experiment B: The concentration of hesperetin was 8.0 µM and the same volume of PPO was added into hesperetin solution 25 times (0→5.0 µM). Two obtained data matrices, D PPO and D hesperetin , were constructed as an extended data matrix [D PPO , D hesperetin ], which was further decomposed using the following equation.
where [C PPO C hesperetin ], S T and [E PPO E hesperetin ] are the matrix for pure component concentration, pure spectra and error, respectively. MCR-ALS analysis steps were described in detail in a previous report [21].

Fluorometric Titration Experiments
Fluorescence spectra were measured at 298, 304 and 310 K by a fluorophotometer (F-7000, Hitachi, Tokyo, Japan). Briefly, the concentration of PPO was 2.0 µM, and the same volume of hesperetin was added into the PPO solution (2.0 mL) 10 times (0→25.0 µM). After each addition of hesperetin, 3 min was timed to equilibrate the mixed system. Then, the emission spectra (290-500 nm) were scanned at fixed excitation wavelength (280 nm) [25]. To deduct the effect of internal filtration, all measured fluorescence data were processed through Formula (5) [26]: where F C and F M represent the corrected and detected fluorescence, respectively. A 1 and A 2 stand for the absorbance values of hesperetin at excitation and emission wavelengths, respectively.
2.2.5. Synchronous Fluorescence Experiments ∆λ, the difference in wavelength between excitation and emission, was set at 15 and 60 nm to investigate the change of peak. The synchronous fluorescence spectra of PPO solution with the addition of hesperetin (0→25.0 µM) were deducted from 240 to 340 nm [27]. The ratios of synchronous fluorescence quenching (RSFQ = 1 − F/F 0 ) were calculated to evaluate the effect of tryptophan (Trp) and tyrosine (Tyr) residues on quenching the fluorescence of PPO [28].

Three-Dimensional (3D) Fluorescence Experiments
The emission spectra of PPO solution with or without hesperetin were recorded once (200-400 nm) when the excitation wavelength changed every 5 nm (200-400 nm). The concentration of PPO and hesperetin in the 2 mL system were 2.0 µM and 15.0 µM, respectively. The obtained data were analyzed and plotted with Matlab software [29].

Molecular Docking
The crystal structure of PPO (PDB ID: 2Y9X) and the 3D molecular structure of hesperetin were downloaded from the RCSB Protein Data Bank and PubChem database, respectively. After the energy minimization of hesperetin and the pretreatment of PPO by removing water, adding hydrogen and adding polarity, the docking of hesperetin and PPO was simulated using Discovery Studio 2016 software (BIOVIA, Beijing, China). Both the top hits and the random conformations were set as 100 times. The best binding conformation owned the lowest interaction energy and was analyzed in details [31].

Molecular Dynamics (MD) Simulation
The more comprehensive information on PPO binding with hesperetin was obtained by MD simulations using GROMACS 5.1 software [32]. Similar to molecular docking, PPO removed water molecules. The topology file of hesperetin went through the pretreatment of the Acpype Server before use. Then, the obtained hesperetin combined with PPO to form the total topology file. The AMBER99SB force field was prepared for hesperetin and PPO. The complex was submerged into a regular dodecahedron water box with a 1 nm edge distance. Na + and Cl − were used for charge neutralization. After the procedures of energy minimization, NVT (300 K), and NPT (1 bar) equilibration, 150 ns MD simulation started to run [33].

Statistical Analysis
Results were expressed as means ± standard deviation. The significant differences (p < 0.05) were considered through one-way ANOVA, followed by Duncan's test using IBM SPSS Statistics 26.

Inhibitory Activity of Hesperetin
To examine the inhibitory ability of hesperetin on monophenolase, L-tyrosine was used as a substrate [10]. As the concentration of hesperetin in the system continuously increased, the corresponding residual activity of PPO was gradually attenuated ( Figure 1A). The IC 50 value of hesperetin was 80.8 ± 1.31 µM, while that of kojic acid was 17.7 ± 0.42 µM, suggesting that the inhibitory ability of hesperetin was weaker than the positive control kojic acid. When the concentration was 100.0 µM, hesperetin decreased the residual activity of PPO to 41.2%, while kojic acid almost completely inhibited the catalytic activity of PPO. Figure 1C,D depict the process curves for L-tyrosine oxidation influenced by different concentrations of hesperetin and kojic acid, respectively. Both hesperetin and kojic acid for reaction rate exhibited a dropping trend with the decreasing concentration ( Figure 1(C-II,1D-II)). As shown in Figure 1(C-III), when hesperetin was at 500.0 µM, the lag time was extended (1.08→10.33 min). Kojic acid at 100.0 µM lengthened the lag time by more than 30 min (Figure 1(D-III)). A previous study reported that the IC 50 value of quercetin on the inhibition of monophenolase activity was 19.4 ± 0.35 µM [10]. Comparatively, hesperetin showed lower inhibitory ability on monophenolase. The structural differences in the compounds may be the important reason for their different inhibitory activities [11]. However, hesperetin was still a valid PPO inhibitor, which may prolong the lag time of PPO oxidation and reduce the reaction rate of the enzyme, consequently inhibiting PPO activity [34].
Regarding previous research [23], L-dopa was prepared as a substrate to determine the inhibitory ability of hesperetin on diphenolase. The inhibition of both hesperetin and kojic acid on PPO increased with their concentrations ( Figure 1B). Once the concentration of hesperetin was greater than 1.0 mM, the residual activity of PPO no longer continued to decline. Similarly, the inhibitory ability of hesperetin on PPO was weaker than kojic acid, whose IC 50 values were 776.0 ± 15.53 µM and 16.2 ± 0.37 µM, respectively. The inhibition of hesperetin on PPO was also weaker than previously studied flavonoids, such as norartocarpetin and luteolin [3]. This may be because hesperetin has less hydroxyl and owns methoxy on its C4 of ring B, with hydrogenation on C2 = C3 [11]. When inhibiting xanthine oxidase's activity, flavonoids also showed a similar regular pattern. Kaempferol had a higher inhibition rate than kaempferide, while their C4 of ring B was converted from hydroxyl to methoxy and the hydrogenation of the C2 = C3 weakened the inhibitory ability of flavonoids; that is, quercetin was superior to taxifolin [28]. Thus, the determination of the inhibitory capacity for hesperetin was in line with previous studies, and further provided theoretical support for developing novel PPO inhibitors.
lines intersected at the origin of the coordinate and showed good linearity. As the concentration of hesperetin increased, the slopes of the lines declined. These results implied that hesperetin reversibly inhibited the activity of PPO, and it did not reduce the amount of PPO available but only attenuated its catalytic rate for L-dopa oxidation [35]. The inhibitory effect of hesperetin against PPO is shown in Figure 1F. Four lines intersected in the second quadrant. Either the slopes or the Y-intercept of four lines augmented with the rising concentrations of hesperetin, which represented the values of Km increased in a concentration-dependent manner while that of Vmax decreased. These re- . c(hesperetin) = 0, 300, 600, 1200 µM for curve a→d, respectively. c(L-dopa) = 0.5 mM; (F) Lineweaver-Burk plots for different concentration of hesperetin (0, 300, 600, 1200 µM for curve a→d, respectively). c(PPO) = 0.2 µM. The secondary plots represent Y-intercept and slope versus [hesperetin], respectively.

Inhibition Kinetics of Hesperetin
At different concentrations of hesperetin, the effect on the enzymatic reaction rate (∆OD/min) by increasing the concentrations of PPO were shown in Figure 1E. The four lines intersected at the origin of the coordinate and showed good linearity. As the concentration of hesperetin increased, the slopes of the lines declined. These results implied that hesperetin reversibly inhibited the activity of PPO, and it did not reduce the amount of PPO available but only attenuated its catalytic rate for L-dopa oxidation [35]. The inhibitory effect of hesperetin against PPO is shown in Figure 1F. Four lines intersected in the second quadrant. Either the slopes or the Y-intercept of four lines augmented with the rising concentrations of hesperetin, which represented the values of K m increased in a concentration-dependent manner while that of V max decreased. These results reflected that hesperetin was the mixed-type inhibitor of PPO [36]. However, a previous study reported that hesperetin was a competitive PPO inhibitor [17]. The inhibition mechanism of hesperetin against PPO deserved further exploration. In addition, the graphs of both Y-intercept and slope versus hesperetin had a good linear relationships, implying that when hesperetin inhibited PPO, there was only one or a class site. According to Equations (2) and (3), the K i and α values were 569.4 ± 6.24 µM and 2.93, respectively. Thus, hesperetin may have a better chance to bind with free PPO instead of PPO-hesperetin complex [37].

Fluorescence Spectra of PPO-Hesperetin Interaction
PPO showed a characteristic fluorescence emission peak at 340 nm. With the gradual addition of hesperetin, the peak red-shifted to 346 nm (Figure 2A). In contrast, the fluorescence emission peak intensity of hesperetin was weak and it was located at 311 nm.
With the continuous addition of PPO, a new fluorescence peak appeared at 340 nm and progressively increased in intensity ( Figure 2B). The position of peak change may be that PPO interacted with hesperetin, and formed a PPO-hesperetin complex. The extended spectral data matrix [D PPO , D hesperetin ] was further analyzed using the MCR-ALS algorithm to judge whether the PPO-hesperetin complex was formed or not [24].

MCR-ALS Decomposition
The decomposition of [D PPO , D hesperetin ] was performed using the singular value decomposition (SVD) model to determine the number of factors associated with chemical species. The first four eigenvalues were 112.75, 42.78, 14.56 and 1.34, implying that three independent factors existed in this complex solution: PPO, hesperetin and PPO-hesperetin complex [24]. As shown in Figure 2C, the resolved profiles of PPO and hesperetin acquired from MCR-ALS analysis (solid lines) overlapped well with the measured profiles (dashed lines), suggesting that the MCR-ALS model correctly resolved the concentration profiles of the components. Either continuously adding hesperetin into the PPO solution ( Figure 2D) or continuously adding PPO into the hesperetin solution ( Figure 2E), the added hesperetin did not react completely with PPO (and vice versa), so that free hesperetin or PPO existed in the system; the concentration profile of the PPO-hesperetin complex always showed an increasing trend. These results provided powerful evidence that PPO interacted with hesperetin and formed PPO-hesperetin complexes [21].

Intrinsic Fluorescence of PPO Quenched by Hesperetin
The binding between PPO and hesperetin could be reflected by the change of fluorescence spectra. PPO had aromatic amino acids (Trp, Tyr, etc.) so that there was a characteristic fluorescence emission peak at 340 nm. Under the same conditions, hesperetin alone only had a slight peak. The fluorescence intensity of PPO showed a regular downward trend with the gradual addition of hesperetin, suggesting that the interaction occurred between hesperetin and PPO. Furthermore, the position of PPO's maximum emission peak conspicuously shifted (340→349 nm), indicating that hesperetin may affect the microenvironment of luminescent amino acid residues in PPO ( Figure 3A) [38].

Mechanism of Fluorescence Quenching
The Stern-Volmer curves for PPO fluorescence quenched by hesperetin at each temperature were obtained through Equation (6) [39]: where F and F 0 denote the fluorescence intensity of PPO with and without hesperetin, respectively. K sv , K q and [Q] represent the fluorescence quenching constant, quenching rate constant and the concentration of hesperetin, respectively. τ 0 represents the average fluorophore lifetime without hesperetin, and its value is usually equal to 10 −8 s [40]. When the temperatures were 298 K, 304 K and 310 K, the K sv values were (1.71 ± 0.01) × 10 4 , (1.76 ± 0.02) × 10 4 and (2.09 ± 0.02) × 10 4 L mol −1 , respectively (Table 1). At the measured temperatures, the order of magnitude of K sv was all 10 4 , meaning that the K q values were all higher than the maximum diffusion collision quenching constant (2.0 × 10 10 L mol −1 s −1 ) [10]. This result signified that PPO bound with hesperetin and formed a static complex, and hesperetin statically quenched PPO's endogenous fluorescence. The K sv values were positively correlated with temperatures, signifying that the increase in temperature favored the fluorescence-quenching effect of hesperetin on PPO [23]. Table 1. Quenching constants K sv , binding constants K a , number of binding sites n and thermodynamic parameters.

Thermodynamic Parameters and Binding Forces
Based on the following equation, the information about the interaction between hesperetin and PPO were obtained, including binding constant (K a ) and binding number (n) [41]: where [Q t ] and [P t ] stand for the concentrations of hesperetin and PPO. At the experimental temperatures, all lines had a good linear relationship. The K a values were 1.98 × 10 4 , 2.05 × 10 4 and 2.22 × 10 4 L mol −1 , which were all in the order of 10 4 L mol −1 , suggesting that when hesperetin bound to PPO, there was a moderate affinity. The higher the temperature, the higher the K a values, indicating that the stability of PPO-hesperetin complex may be disturbed by rising temperature [23]. The n values at 298, 304 and 310 K were all equal to 1, implying that only one site existed when hesperetin interacted with PPO. The values of ∆H • , ∆S • and ∆G • , which refer to enthalpy change, entropy change and free energy change, respectively, were obtained by the van 't Hoff equations [42]: where T represents temperature; R represents gas constant, and its value is 8.314 J mol −1 K −1 . The negative ∆G • values indicated that the combining action of hesperetin and PPO was spontaneous [21]. The positive ∆H • value (17.47 KJ mol −1 ) meant that the binding of hesperetin and PPO was an endothermic process. The ∆S • value (107.22 J mol −1 K −1 ) was also positive ( Table 1), indicating that the interaction between hesperetin and PPO made the PPO-hesperetin system become more disordered. A previous study in which hesperetin interacted with glutenin reported the same results [43]. Furthermore, based on the well-known theory of Ross and Subramanian [44], the positive ∆H • and ∆S • values show that hesperetin was combined with PPO mainly by hydrophobic interaction [33].

Conformational Changes of PPO
The synchronous fluorescence spectra were available for analyzing how hesperetin affected the microenvironment of Tyr and Trp residues in PPO ( Figure 3B,C). Their fluorescence intensities were consistently attenuated by the gradual addition of hesperetin (0→25.0 µM). The position of the maximum emission peak for Trp residues had a slight shift (285→287 nm), indicating that hesperetin augmented the polarity of the microenvironment around Trp residues in PPO, while it had no effect on that around Tyr residues. Furthermore, Tyr residues had higher values of RSFQ compared with Trp residues (Figure 3D), implying that Tyr residues may be closer to their binding site, and contribute more when hesperetin quenched PPO's fluorescence [45].
Another method to judge whether hesperetin changed the conformation of PPO was to compare the 3D fluorescence spectra of PPO and PPO-hesperetin ( Figure 3E,F). Peak a (λ ex = λ em ) and Peak b (2 λ ex = λ em ) denote the Rayleigh scattering peak and second-order scattering peaks, respectively. Peak 1, the feature of Trp and Tyr residues, whose intensity was weakened (601.2→296.9) and position had a movement (λ ex /λ em = 280/340 nm→275/350 nm). The position of peak 2 (the characteristic peak of polypeptide backbone structure generated through π→π* transition) shifted from 230/335 nm to 235/340 nm, and its intensity dropped from 383.8 to 65.9. The alteration of peak 1 further implied that hesperetin influenced Trp and Tyr residues. The change of Peak 2 meant that the addition of hesperetin may change the conformation of PPO as well [46].
The CD spectrum is so sensitive to the change of secondary structure of the protein that is a powerful means to further detect the influence of hesperetin on PPO's conformation. As shown in Figure 3G, the two negative absorption bands at 208 nm and 222 nm were considered as the characteristic CD spectrum of α-helix structure [47], whose content was positively correlated with its intensity. Thus, hesperetin increased the content of α-helix in PPO [30]. As shown in Figure 3H, with the increase in [hesperetin] versus [PPO] (0:1→1:1→4:1), the contents of α-helix (33.37%→35.19%→39.68%) and β-turn (21.12%→22.46%→24.40%) increased, while the content of β-sheet reduced (21.12%→17.95%→13.39%). This phenomenon meant that hesperetin might induce the structure of PPO to partially contract, and the tighter conformation of the enzyme may make the entry of substrates into the active center more difficult, thus inhibiting PPO activity [18]. In addition to chelating the copper ions in the active center to inhibit PPO activity [17], hesperetin changed PPO's conformation, which may also be one of the reasons for the inhibition of hesperetin against PPO.

Molecular Docking
The information about the binding site, binding mode and force type for the interaction between hesperetin and PPO could be intuitively observed by molecular docking [48]. Hesperetin entered the cavity of PPO and interacted with surrounding hydrophobic amino acid residues ( Figure 4A). Hesperetin bound near the active center (containing two Cu 2+ ) of PPO ( Figure 4B). The A7-OH of hesperetin formed a hydrogen bond with Met280, whose distance was 4.64 Å. The other two hydrogen bonds were formed between the carbonyl in the C ring of hesperetin and His89 and His259 and the distances were 4.39 Å and 5.79 Å, respectively ( Figure 4C,D). Moreover, hesperetin interacted with Val283, Phe264, His85, Asn260, Val248 and His263 via hydrophobic forces, while with Phe292, His61, Phe90, Glu256, His244, Asn260, Phe264 and Gly281 via van der Waals forces. Thus, hydrogen bonds, hydrophobic interactions and van der Waals forces were all important forces in the formation of the PPO-hesperetin complex. Si et al. [17] also simulated the molecular docking of hesperetin with PPO using two different docking programs. They found that Met280, His61, His85 and His259 in PPO were important sites in forming hydrogen bonds with hesperetin. The findings in our study were similar to the report. It is noteworthy that hesperetin interacted with two copper ions and the critical residues around the active center of PPO. This phenomenon was also found by Song et al. [23], who reported that the interactions between the compound and PPO may influence the stabilization of the dinuclear copper active center of PPO.
All in all, hesperetin entered PPO's hydrophobic cavity, bound near its active center, acted as an obstruction to occupy the substrate's binding site in PPO, prevented the substrate being catalyzed by PPO and may change the structure of PPO as well, thus leading to decreased PPO activity [49].

Molecular Dynamics
Experiments can only give macroscopic results, while MD simulations can provide atomic-level insights into the binding of hesperetin to PPO. Therefore, this method can be used to supplement and verify the known results.
The root mean square deviation (RMSD) value represents the stability of the systems [23]. By comparing free PPO with PPO-hesperetin complex in RMSD values, we could determine whether hesperetin affects the stability of PPO. The more stable systems have less variation in RMSD values. The RMSD values of PPO and PPO-hesperetin systems tended to increase, reached a relatively stable state after 100 ns, and finally stabilized at around 0.22 nm ( Figure 5A). The RMSD values of PPO-hesperetin had less volatility in the first 30 ns, which may be caused by the formation of hydrogen bonds in the complex system [50]. After stabilization, the complex owned the lower RMSD values, implying that adding hesperetin into the system improved its stability [51]. The radius of gyration (Rg) value is negatively correlated with the structural denseness of the enzyme [52]. The Rg value of the PPO-hesperetin system remained at around 2.065 nm, which was smaller than the free PPO system at 2.075 nm ( Figure 5B), suggesting that hesperetin may strengthen the tightness of the PPO structure [42].
The root mean square fluctuation (RMSF) value reflects the flexibility of amino acid residues in PPO influenced by hesperetin [53]. Overall, the free PPO system had a higher  The radius of gyration (Rg) value is negatively correlated with the structural denseness of the enzyme [52]. The Rg value of the PPO-hesperetin system remained at around 2.065 nm, which was smaller than the free PPO system at 2.075 nm ( Figure 5B), suggesting that hesperetin may strengthen the tightness of the PPO structure [42].
The root mean square fluctuation (RMSF) value reflects the flexibility of amino acid residues in PPO influenced by hesperetin [53]. Overall, the free PPO system had a higher RMSF value than the complex system; in particular, the residues 240-260 ( Figure 5C). This result suggested that hesperetin could improve the stability of the PPO structure by restricting some amino acids, which may be involved in the binding process of hesperetin with PPO.
The solvent-accessible surface area (SASA) value reflects the change in the surface area between the system and the solvent, and its value increases with the decreasing hydrophobicity of the enzyme [23]. The decreased SASA value ( Figure 5D) indicated that hesperetin reduced PPO's hydrophilicity and tightened its structure [54]. Hesperetin only had a feeble effect on the microenvironment of Trp and Tyr in PPO ( Figure 5E,F). All the MD simulation results were consistent with the results from fluorescence and CD spectra experiments.

Conclusions
Hesperetin could inhibit the oxidation of L-tyrosine and L-dopa catalyzed by PPO, which was a reversible and mixed-type PPO inhibitor. MCR-ALS analysis showed that PPO interacted with hesperetin and formed the PPO-hesperetin complex. The combination of hesperetin and PPO was mainly driven by hydrophobic interactions, while only one site existed for their binding. The addition of hesperetin affected the microenvironment of residues in PPO and increased PPO's α-helix content, which represented the tighter molecular structure of PPO, and thus may make the entry of the substrate into, and the release of the catalytic products from, the active center more complex, leading to decreased PPO activity. This study may provide novel insights into the inhibition mechanism of hesperetin against PPO and theoretical guidance for the development of new PPO inhibitors. However, the anti-browning effect of hesperetin in real food systems need be further investigated.

Data Availability Statement:
The data are available from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.

PPO
polyphenol oxidase IC 50 half-maximum inhibitory concentration MCR-ALS multivariate curve resolution-alternate least squares CD circular dichroism SVD singular value decomposition RSFQ ratios of synchronous fluorescence quenching 3D three-dimensional MD molecular dynamics RMSD root mean square deviation Rg radius of gyration RMSF root mean square fluctuation SASA solvent accessible surface area